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Active adaptive selective control system 





(57) An active adaptive control system introduces a 
control signal from an output transducer (1 4) to combine 
with the system input signal (6) and yield a system out- 
put signal (8). An error transducer (16) senses the sys- 
tem output signal and provides an error signal (44). An 
adaptive filter model (40) has a model input from a ref- 
erence signal (42) correlated to the system input signal, 
and an output outputting a correction signal (46) to the 



output transducer to introduce the control signal. Per- 
formance of the model is selectively controlled to control 
the signal sent to the output transducer. Various moni- 
toring and control methods are provided, including spec- 
tral leak signal monitoring and control, correction signal 
monitoring and control, frequency responsive spectral 
transfer function processing of the leak signal and/or the 
correction signal, reference signal processing, and 
fuzzy logic control. 
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Description 

BACKGROUND AND SUMMARY 

The invention relates to active adaptive control sys- 
tems, and more particularly to improvements for selec- 
tively controlling performance of the active adaptive 
model. 

The invention arose during continuing development 
efforts relating to the subject matter of U.S. Patents 
4,837,834, 5,172,416, 5,278,913, 5,386,477, 
5,390,255, 5,396,561, and co-pending U.S. Applica- 
tions Serial No. 08/166,698, filed December 14, 1993, 
(EP-A-0 661 807), 08/247,561 , filed May 23, 1 994 (EP- 
A-0 684 594), 08/264 ; 510 ; filed June 23, 1994, 
08/340,613, filed November 16, 1994, 08/369,925, filed 
January 6, 1995, (EP-A-0 721 179) incorporated herein 
by reference. 

Active acoustic attenuation involves injecting a can- 
celing acoustic wave to destructively interfere with and 
cancel an input acoustic wave. In an active acoustic at- 
tenuation system, the output acoustic wave is sensed 
with an error transducer, such as a microphone or an 
accelerometer, which supplies an error signal to an 
adaptive filter control model which in turn supplies a cor- 
rection signal to a canceling output transducer or actu- 
ator, such as a loudspeaker or a shaker, which injects 
an acoustic wave to destructively interfere with the input 
acoustic wave and cancel or reduce same such that the 
output acoustic wave at the error transducer is zero or 
some other desired value. 

An active adaptive control system minimizes an er- 
ror signal by introducing a control signal from an output 
transducer to combine with the system input signal and 
yield a system output signal. The system output signal 
is sensed with an error transducer providing the error 
signal. An adaptive filter model has a model input from 
a reference signal correlated with the system input sig- 
nal, an error input from the error signal, and outputs a 
correction signal to the output transducer to introduce a 
control signal matching the system input signal, to min- 
imize the error signal. The filter coefficients are updated 
according to a weight update signal which is the product 
of the reference signal and the error signal. 

The present invention is applicable to active adap- 
tive control systems, including active acoustic attenua- 
tion systems. In one embodiment, the invention maxi- 
mizes model performance by concentrating model ad- 
aptation in frequency ranges of interest and protects the 
output transducer against overdriving of same. Perform- 
ance of the model is spectrally controlled to maximize 
the correction signal sent to the output transducer such 
that at frequencies where maximum power from the out- 
put transducer reaches the error transducer, the correc- 
tion signal supplied to the output transducer is maxi- 
mized, and at frequencies where minimum power from 
the output transducer reaches the error transducer, the 
correction signal supplied to the output transducer is 



minimized. This maximizes model performance by con- 
centrating model adaptation on the portion of the input 
signal which the model can control or cancel or where 
it is desired to do so, and constrains model adaptation 
5 as to those portions of the input signal which it cannot 
cancel or control or where it is not desired to do so. The 
latter is desired for stability of the model algorithm where 
active control solutions sometimes require more actua- 
tor power than is available or desirable. Actuators, am- 
plifiers, etc. have limitations that adversely affect control 
algorithms. Pushed beyond capacity, the control output 
or power available from the secondary source or output 
transducer may exhibit saturation, clipping, or otherwise 
nonlinear behavior. Excessive control effort can result 
in damaged actuators, excessive power consumption, 
and instability in the control algorithm. 

It is known in the prior art to provide weight update 
signal leakage to counteract the adaptive process. This 
is done by implementing an exponential decay of the 
filter coefficients, intentionally defeating control effort, 
Widrowand Stearns, Adaptive Signal Processing , Pren- 
tice-Hall, Inc., Engelwood Cliffs, NJ, 1985, pages 
376-378. The exponential decay is typically selected to 
be slow such that the adaptive process toward a control 
solution dominates. A deficiency of this method is that 
it unilaterally, across all power levels and frequencies, 
degrades performance. Such leakage is useful for lim- 
iting control effort and enhancing numerical stability but 
performance suffers because of the lack of considera- 
tion for regions where the control effort is in an accept- 
able range. 

In other embodiments of the invention, model ad- 
aptation is selectively controlled to provide desired per- 
formance. Model performance is controlled by fuzzy log- 
ic to provide self-designing control architecture using 
fuzzy rules and/or to control a filter transfer function and/ 
or to control filter weights used in an update process for 
feedforward and/or feedback, including FIR (finite im- 
pulse response) and MR (infinite impulse response) ap- 
plications and/or to control magnitude and/or rate of 
change of a leak signal degrading performance of the 
model. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a schematic illustration of an active adap- 
tive control system known in the prior art. 

Fig. 2 is a schematic illustration of an active adap- 
tive control system in accordance with co-pending U.S. 
Application S.N. 08/264,510, filed June 23, 1994. 

Fig. 3 is a graph showing performance of the system 
of Fig. 2. 

Fig. 4 is a graph further showing performance of the 
system of Fig. 2. 

Fig. 5 is a graph showing an alternate performance 
of the system of Fig. 2. 

Fig. 6 is a graph further showing alternate perform- 
ance of the system of Fig. 2. 
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Fig. 7 is a schematic illustration of an active adap- 
tive control system in accordance with the present in- 
vention. 

Fig. 8 is a graph showing frequency versus output 
and illustrates performance of the system of Fig. 7. 

Fig. 9 is a schematic illustration of an active adap- 
tive control system. 

Fig. 1 0 is a schematic illustration of an adaptive filter 
model, and illustrates a principle employed by a system 
in accordance with copending U.S. Application Serial 
No. 08/166,698, filed December 14, 1993. 

Fig. 1 1 is like Fig. 10 and shows another manner of 
implementing the principle thereof. 

Fig. 1 2 is a schematic illustration of an active adap- 
tive control system in accordance with the system of the 
'698 application. 

Fig. 13 shows a further embodiment of the system 
of Fig. 12. 

Fig. 14 shows a further embodiment of the system 
of Fig. 12. 

Fig. 1 5 is a schematic illustration of an active adap- 
tive control system in accordance with the present in- 
vention. 

Fig. 16 shows a further embodiment of the system 
of Fig. 15. 

Fig. 17 shows a further embodiment of the system 
of Fig. 15. 

Fig. 1 8 is a schematic illustration of an active adap- 
tive control system in accordance with the present in- 
vention. 

Fig. 19 shows a further embodiment of the system 
of Fig. 18. 

Fig. 20 shows a further embodiment of the system 
of Fig. 18. 

Fig. 21 shows a further embodiment of the system 
of Fig. 18. 

Figs. 22-24 illustrate a further embodiment of the 
invention. 

DETAILED DESCRIPTION 

Fig. 1 shows an active adaptive control system sim- 
ilar to that shown in U.S. Patent 4,677,676, incorporated 
herein by reference, and uses like reference numerals 
therefrom where appropriate to facilitate understanding. 
The system introduces a control signal from a second- 
ary source or output transducer 14, such as a loud- 
speaker shaker, or other actuator or controller, to com- 
bine with the system input signal 6 and yield a system 
output signal 8. An input transducer 10, such as a mi- 
crophone, accelerometer, or other sensor, senses the 
system input signal and provides a reference signal 42. 
An error transducer 16, such as a microphone, acceler- 
ometer, or other sensor, senses the system output sig- 
nal and provides an error signal 44. Adaptive filter model 
40 adaptively models the system and has a model input 
from reference signal 42 correlated to system input sig- 
nal 6, and an output outputting a correction signal 46 to 



output transducer 14 to introduce the control signal ac- 
cording to a weight update signal 74. Reference signal 
42 and error signal 44 are combined at multiplier 72 to 
provide the weight update signal through delay element 

s 73. In a known alternative, the reference signal 42 may 
be provided by one or more error signals, in the case of 
a periodic system input signal, "Active Adaptive Sound 
Control In A Duct: A Computer Simulation", J.C. Bur- 
gess, Journal of Acoustic Society of America, 70(3), 

10 September 1981, pages 715-726, U.S. Patents 
5,206,911 , 5,216,722, incorporated herein by reference. 

In updating the filter coefficients, and as is standard, 
one or more previous weights are added to the current 
product of reference signal 42 and error signal 44 at 

15 summer 75. As noted above, it is known in the prior art 
to provide exponential decay of all of the filter coeffi- 
cients in the system. Leakage factor y at 77 multiplies 
one or more previous weights, after passage through 
one or more delay elements 73, by an exponential decay 

20 factor less than one before adding same at summer 75 
to the current product of reference signal 42 and error 
signal 44, Adaptive Signal Processing , Widrow and 
Stearns, Prentice-Hall, Inc., Engelwood Cliffs, NJ, 1985, 
pages 376-378, including equations 13.27 and 13.31. 

25 As noted above, a deficiency of this method is that it 
reduces control effort and degrades performance 
across all power levels, regardless of whether such re- 
duced effort is desired. 

In the '510 application, selective leakage of the 

30 weight update signal is provided in response to a given 
condition of a given parameter, to control performance 
of the model on an as needed basis. In the preferred 
embodiment, leakage is varied as a function of correc- 
tion signal 46. A variable leakage factory is provided at 

35 79 in Fig. 2, replacing fixed y 77 of Fig. 1 . Leakage factor 
yat 79 is varied from a maximum value of 1 .0 affording 
maximum control effort, to a minimum value such as ze- 
ro providing minimum control effort. 

Leakage is varied as a function of the output power 

40 of correction signal 46 supplied from the output of model 
40 to output transducer 1 4. In the embodiment in Fig. 3, 
the leakage is varied as a discontinuous step function 
of the output power of the correction signal. When the 
output power exceeds a given threshold at 81, y is 

45 abruptly, nonlinearly changed as a step function from a 
first level 83 to a second level 85. The reduction at 85 
reduces the weight update signal summed at summer 
75 with the product of the reference signal 42 and error 
signal 44 from multiplier 72, and hence reduces the 

50 weight update signal supplied to model 40. The noted 
reduction of yat threshold 81 increases leakage of the 
weight update signal, Fig. 4, from level 87 to level 89. 

In another embodiment as shown in Fig. 5, leakage 
is varied as a continuous function of the output power 

55 of the correction signal. In Fig. 5, y is maintained at level 
83 until output power reaches threshold 81 , and then is 
linearly decreased as shown at 91 as a continuous lin- 
early changing value as a function of increasing output 



3 



5 



EP 0 759 606 A2 



6 



power above threshold 81 . As shown in Fig. 6, leakage 
is maintained at level 87 until output power reaches 
threshold 81, and then is linearly increased at 93 as a 
continuous linearly changing value as a function of in- 
creasing output power above threshold 81 . 

Other variations of leakage are possible for provid- 
ing selective leakage of the weight update signal to de- 
grade performance of the model. The leakage is adjust- 
ably varied to vary performance of the model by multi- 
plying a previous weight update value by variable y 79 
and adding the result at summer 75 to the product of 
reference signal 42 and error signal 44 from multiplier 
72. y79 is varied as a function of correction signal 46, 
preferably the output power of such correction signal. 

Fig. 7 illustrates the present invention and uses like 
reference numerals from Fig. 2 and from Figs. 19 and 
20 of the incorporated '676 patent. The transfer function 
from output transducer 1 4 to error transducer 1 6 is mod- 
eled with an adaptive filter model C at 142 ; as in the 
incorporated '676 patent. Auxiliary random signal 
source 140 introduces a random signal into the output 
of model 40 at summer 1 52 and into the C model at 1 48. 
The auxiliary random signal from source 140 is random 
and uncorrected with the system input signal 6 and in 
preferred form is provided by a Galois sequence, M.R. 
Schroeder, "Number Theory In Science And Communi- 
cations", Berlin, Springer-Berlag, 1984, pages 252-261 , 
though other random uncorrelated signal sources may 
be used. The Galois sequence is a pseudo random se- 
quence that repeats after 2 M -1 points, where M is the 
number of stages in a shift register. The Galois se- 
quence is preferred because it is easy to calculate and 
can easily have a period much longer than the response 
time of the system. The input 1 48 to C model 1 42 is mul- 
tiplied with the error signal from error transducer 16 at 
multiplier 68, and the resultant product provided as 
weight update signal 67. Model 142 models the transfer 
function from output transducer 14 to error transducer 
16, including the transfer function of each. Alternatively, 
the transfer function from output transducer 14 to error 
transducer 1 6 may be modeled without a random signal 
source, as in U.S. Patent 4,987,598, incorporated here- 
in by reference. Auxiliary source 140 introduces an aux- 
iliary random signal such that error transducer 16 also 
senses the auxiliary signal from the auxiliary source. 
The auxiliary signal may be introduced into the recursive 
loop of the A and B filters as in Fig. 1 9 of the incorporated 
'676 patent at summer 1 52, or alternatively the auxiliary 
signal may be introduced into the model after the recur- 
sive loop : i.e. introducing the auxiliary signal only to line 
46, and not to line 47. A copy of model 1 42 is provided 
at 1 45 to compensate the noted transfer function, as in 
the incorporated '676 patent. 

As in the '510 application, the present system and 
method involves introducing a control signal from output 
transducer 1 4 to combine with system input signal 6 and 
yield system output signal 8, sensing the system output 
signal with error transducer 16 and providing an error 



signal 44, providing adaptive filter model 40 having a 
model input from reference signal 42 correlated to sys- 
tem input signal 6, and an output outputting a correction 
signal 46 to output transducer 14 to introduce the control 
5 signal according to weight update signal 74. In the 
present invention, performance of model 40 is spectrally 
controlled to maximize the signal sent to output trans- 
ducer 14 at frequencies of interest or where it can be 
maximally effective, and minimize the signal sent to out- 
put transducer 1 4 at frequencies of non interest or where 
it is only minimally effective or is ineffective. A spectral 
leak signal is provided which degrades performance of 
the model. The leak signal is controlled according to fre- 
quency. In Fig. 7, the correction signal from the output 
of model 40 is monitored, and the leak signal is control- 
led in response thereto. The correction signal is filtered 
by filter 95 to provide the leak signal. The correction sig- 
nal from the output of model 40 is spectrally processed 
by supplying the correction signal through the frequency 
responsive spectral transfer function provided by filter 
95 to provide the leak signal to the error input of model 
40. In Fig. 7, the correction signal is spectrally proc- 
essed such that at frequencies where maximum power 
from output transducer 14 reaches error transducer 16, 
the correction signal is maximized, and at frequencies 
where minimum powerfrom output transducer 14 reach- 
es error transducer 16, the correction signal is mini- 
mized. 

The transfer function between output transducer 1 4 
and error transducer 16 is modeled with C model 142. 
Correction signal 46 is spectrally processed by a func- 
tion of model 142. In Fig. 7, such function is the inverse 
of C model 142 as provided at inverse C model, C' 1 , at 
95. The output of inverse C model 95 is supplied to an 
optional peak detector as provided by summer 97 com- 
paring the output of inverse C model 95 with a desired 
peak value 99. When the output of model 95 rises above 
level 99, the positive output of summer 97 controls var- 
iable leakage factor y at 79, as above. The inverse C 
model includes the inverse of the transfer function of 
output transducer 14, inverse S, S" 1 , and/or the inverse 
of the transfer function of the error path, inverse E, E _1 , 
between output transducer 14 and error transducer 16. 

In addition to or in place of peak detector 97, control 
logic may be used to respond to the output of inverse C 
model 95 and control leakage factor y at 79 according 
to designated conditions or rules to generate or compute 
a leak value or control the leaking process ; to be further 
described hereinafter. In other embodiments, filter 95 is 
a displacement function of output transducer 1 4 such as 
a loudspeaker, to protect the latter against overdriving. 
In other embodiments, an RMS (root mean square) to 
DC (direct current) conversion function is provided be- 
tween filter 95 and peak detector 97 for average level 
detection to control convergence rate of the leak proc- 
ess. The transfer function of filter 95 may be linear or 
nonlinear. 

As seen in Fig. 8, at frequencies where maximum 
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power from output transducer 14 reaches error trans- 
ducer 16, as at frequency regions 101 and 103, C model 
142 has a maximum transfer characteristic, and the in- 
verse of the C model at 95 has a minimum transfer char- 
acteristic, as shown at 1 05 and 1 07. The minimum trans- 
fer characteristic at 105 and 107 minimizes leakage of 
the update signal to model 40, to enable maximum out- 
put of model 40. At frequencies where minimum power 
from output transducer 14 reaches error transducer 16, 
as at region 109, C model 142 has a minimum transfer 
characteristic, and the inverse of the C model at 95 has 
a maximum transfer characteristic, as shown at 1 1 1 . The 
maximum transfer characteristic at 1 1 1 maximizes leak- 
age of the update signal, to minimize the output of model 
40. Inverse C model 95 spectrally senses correction sig- 
nal 46 and provides selective leakage of weight update 
signal 74 in response thereto, to control performance of 
model 40 according to frequency, to optimize perform- 
ance of model 40 in frequency ranges such as 1 01 and 
103 where the model can effectively control or cancel 
input signal 6. Outside of such ranges, e.g. at 218, in- 
verse C model 95 minimizes performance of model 40, 
to avoid using computation or adaptation power where 
inefficient or unneeded ; and to prevent the model from 
continually trying to generate an output in regions where 
it is ineffective to attempt to control or cancel input signal 
6. Model 95 affords frequency weighting of the weight 
update signal. 

Fig. 9 is similar to Fig. 5 of U.S. Patent 4,677,676, 
incorporated herein by reference, and uses like refer- 
ence numerals to facilitate understanding. The system 
introduces a control signal from an output transducer 1 4, 
such as a loudspeaker shaker, or other actuator or con- 
troller, to combine with the system input signal 6 and 
yield a system output signal 8. An input transducer 10, 
such as a microphone, accelerometer, or other sensor, 
senses the system input signal. An error transducer 1 6, 
such as a microphone, accelerometer, or other sensor, 
senses the system output signal and provides an error 
signal. Adaptive filter model 40 adaptively models the 
system and has a model input 42 from input transducer 
10, an error input 44 from error transducer 16, and a 
model output 46 outputting a correction signal to output 
transducer 1 4 to introduce the control signal. In a known 
alternative, the input signal at 42 may be provided by 
one or more error signals, in the case of a periodic sys- 
tem input signal, "Active Adaptive Sound Control In A 
Duct: A Computer Simulation", J.C. Burgess, Journal of 
Acoustic Society of America, 70(3), September 1981, 
pages 715-726, U.S. Patents 5 ; 206 ; 911 , 5,216,722, in- 
corporated herein by reference. 

The system of the '698 application provides an ac- 
tive adaptive control system wherein the performance 
of model 40 is intentionally and selectively constrained 
by driving the output 46 of the model towards zero in 
response to a given condition of a given parameter. For 
example, in an active noise control system, it may be 
desirable to cancel noise only in a given frequency band, 



and leave the noise uncanceled for frequencies outside 
the band. In other control applications, it may be desir- 
able to selectively control the system output signal by 
selectively controlling introduction of the control signal 
5 from output transducer 1 4 to match or not match the sys- 
tem input signal. 

One manner of constraining system performance is 
to drive the output of model 40 towards zero and away 
from a value matching system input signal 6. One way 
of accomplishing this is shown in Fig. 10 wherein the 
output of model 40 is supplied to its error input, such 
that when the model adapts to drive its error input to- 
wards zero, the output of the model is necessarily also 
driven towards zero. Fig. 11 shows another manner of 
implementing this principle wherein a copy of the model 
is provided at 200, and the output of model copy 200 
supplies the error signal to the error input of model 40. 
In Fig. 11, model 40 adapts to drive its error input to- 
wards zero, which in turn requires that the output of copy 
200 be driven towards zero, which in turn means that 
the output of model 40 is driven towards zero because 
M copy 200 is a duplicate of model 40. These principles 
are utilized in the present invention. 

Model 40, Fig. 1 2, normally adapts to a converged 
condition wherein its output at 46 provides a correction 
signal to output transducer 14 which outputs a control 
signal matching the system input signal or a designated 
relative value correlated thereto. For example, in a noise 
cancellation system, the matching control signal from 
output transducer 14 cancels the input noise. In the sys- 
tem of the '698 application, in response to a given con- 
dition of a given parameter, the output of model 40 is 
driven towards zero by driving the output of the model 
towards its error input such that when the model adapts 
to drive the error signal towards zero, the output of the 
model is also driven towards zero. This is accomplished 
by providing a copy 200 of model 40, and supplying the 
output of the copy to an error input 202 of the model 
which is summed at summer 204 with the error signal 
at error input 44 from error transducer 16. The model 
adapts to drive the error input towards zero which in turn 
requires that the output of copy 200 and hence the out- 
put of model 40 are driven towards zero, to provide the 
noted constrained performance. The driving of model 
output 46 towards zero provides a zero or at least a re- 
duced correction signal to output transducer 1 4 to con- 
strain or reduce modification and/or cancellation of the 
system input signal 6. 

A random signal is provided at 206 from an auxiliary 
random signal source 208, preferably provided by a Ga- 
lois sequence, M.R. Schroeder, "Number Theory In Sci- 
ence And Communications", Berlin, Springer-Berlag, 
1984, pages 252-261, though other random signal 
sources may be used, uncorrelated with the system in- 
put signal 6. The Galois sequence is a pseudo random 
sequence that repeats after 2 M -1 points, where M is the 
number of stages in a shift register. The Galois se- 
quence is preferred because it is easy to calculate and 
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can easily have a period much longer than the response 
time of the system. The random signal is supplied 
through a stopband filter 21 0 to model copy 200 at 21 2. 
Stopband filter 21 0 blocks frequencies in the stopband, 
and passes frequencies outside the stopband. This ar- 
rangement provides a spectral leak signal at 202 in re- 
sponse to a given condition of a given parameter, for 
example a frequency outside the stopband. In such im- 
plementation, the noted given parameter is frequency, 
and the given condition is a designated sub-optimum 
performance band outside the stopband. 

The spectral leak error signal at 202 drives the cor- 
rection signal at model output 46 towards zero and pro- 
vides sub-optimum performance of model 40. Outside 
of the sub-optimum performance band, i.e. within the 
stopband of filter 210, there is no signal at 212 and 
hence the output of copy 200 is undefined, and the error 
signal from error transducer 1 6 at error input 44 is max- 
imally effective and model 40 optimally responds thereto 
and drives the correction signal at output 46 toward a 
value matching the system input signal 6. When the 
spectral leak signal is present at error input 202, it con- 
strains performance of model 40 by driving or at least 
attempting to drive the correction signal at model output 
46 towards zero. The relative influence or amplitudes of 
the error signals at error inputs 44 and 202 are adjusted 
to provide the desired relative dominance. Where it is 
desired to eliminate all modification and/or cancellation 
of the system input signal when the frequency is outside 
the stopband of filter 210, then the noted relative ampli- 
tudes are set such that the error signal at error input 202 
dominates the error signal at error input 44, and hence 
the correction signal at model output 46 is driven to- 
wards zero and away from a value matching the system 
input signal 6. 

The method of the '698 application involves driving 
error input 44 to drive the correction signal at model out- 
put 46 toward a value matching the system input signal, 
and selectively driving error input 202 to drive the cor- 
rection signal at model output 46 away from the match- 
ing value by driving the correction signal towards zero. 
The arrangement provides a spectral leaksignal to error 
input 202 in response to the noted given condition of the 
given parameter, e.g. a frequency outside the stopband, 
such that in the presence of the given condition, the 
spectral leak signal drives the correction signal at model 
output 46 towards zero, and in the absence of the given 
condition the error signal at error input 44 drives the cor- 
rection signal at model output 46 towards a value match- 
ing the system input signal 6. 

Stopband filter 210 blocks frequencies in a given 
stopband at which modification or cancellation of the 
system input signal 6 by model 40 is desired. Filter 210 
passes frequencies in a given passband at which mod- 
ification or cancellation of the system input signal by 
model 40 is undesired or not possible. The control signal 
output by output transducer 1 4 is driven toward a value 
matching the system input signal 6 only for frequencies 



in the stopband. At frequencies in the stopband, the er- 
ror signal at error input 44 is dominant, and the control 
signal output by output transducer 14 is driven toward 
a value matching the system input signal 6. At frequen- 
5 cies in the passband, the error signal at error input 202 
is dominant, and the control signal output by output 
transducer 14 is driven away from a value matching the 
system input signal 6. 

Fig. 1 3 is similar to Figs. 1 9 and 20 of the incorpo- 
rated '676 patent, and uses like reference numerals 
where appropriate to facilitate understanding. As noted 
in the incorporated '676 patent, model M at 40 is pref- 
erably an adaptive recursive filter having a transfer func- 
tion with both poles and zeros. Model M is provided by 
an MR, infinite impulse response, filter, e.g. a recursive 
least mean square, RLMS, filter having a first algorithm 
filter provided by an FIR, finite impulse response, filter, 
e.g. a least mean square, LMS, filter A at 12, and a sec- 
ond algorithm filter provided by an FIR fi Iter, e.g. an LMS 
filter, B at 22. Filter A provides a direct transfer function, 
and filter B provides a recursive transfer function. The 
transfer function from output transducer 14 to error 
transducer 16 is modeled by a filter e.g. an LMS or 
RLMS filter, C at 1 42, as in the incorporated '676 patent. 

Auxiliary random signal source 140 introduces a 
random signal into the output of model 40 at summer 
152 and into the C model at 148. The auxiliary random 
signal from source 140 is random and uncorrelated with 
the system input signal 6, and is also uncorrelated with 
auxiliary random signal source 208, and in preferred 
form is provided by a Galois sequence, M.R. Schroeder, 
"Number Theory In Science And Communications", Ber- 
lin, Springer-Berlag, 1 984, pages 252-261 , though other 
random uncorrelated signal sources may be used. The 
Galois sequence is a pseudo random sequence that re- 
peats after 2 M -1 points, where M is the number of stages 
in a shift register. The Galois sequence is preferred be- 
cause it is easy to calculate and can easily have a period 
much longer than the response time of the system. The 
input 1 48 to C model 1 42 is multiplied with the error sig- 
nal from error transducer 1 6 at multiplier 68, and the re- 
sultant product provided as weight update signal 67. 
Model 142 models the transfer function from output 
transducer 14to error transducer 16, including the trans- 
fer function of each. Alternatively, the transfer function 
from output transducer 14 to error transducer 16 may 
be modeled without a random signal source, as in U.S. 
Patent 4,987,598, incorporated herein by reference. 
Auxiliary source 1 40 introduces an auxiliary random sig- 
nal such that error transducer 1 6 also senses the auxil- 
iary signal from the auxiliary source. The auxiliary signal 
may be introduced into the recursive loop of the A and 
B filters as in Fig. 19 of the incorporated '676 patent at 
summer 1 52 ; or alternatively the auxiliary signal may be 
introduced into the model after the recursive loop, i.e. 
introducing the auxiliary signal only to line 46, and not 
to line 47. As in the incorporated '676 patent, copies of 
model 142 are provided at 144 and 146 to compensate 
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the noted transfer function. 

The outputs of filters A and B are summed at sum- 
mer 48, whose output is summed at summer 152 with 
the output of random signal source 140 to provide an 
output resultant sum which provides the model output 
at 46 supplying the noted correction signal to output 
transducer 14. The output of model 142 is summed at 
summer 64 with the output of error transducer 16, and 
the resultant sum supplied as the error input to model 
1 42 and as an error input to model 40. Alternatively the 
output of error transducer 16 may be supplied directly 
to an error input of model 40 without being supplied 
through summer 64. 

M copy 200, Fig. 1 2, is provided by a copy of A filter 
12 at A copy 214, Fig. 13, and a copy of B filter 22 at B 
copy 216. Multiplier 218 multiplies the output 220 of A 
copy 21 4 and the input 222 of A copy 214, and supplies 
the output resultant product at 224 to summer 226. Mul- 
tiplier 228 multiplies the output 230 of C copy 144 and 
the error input at 44, and supplies the output resultant 
product at 232 to summer 226. Summer 226 sums the 
inputs 232 and 224, and supplies the output resultant 
sum as weight update signal 74 to A filter 12. Multiplier 
234 multiplies the output 236 of B copy 21 6 and the input 
238 of B copy 216, and supplies the output resultant 
product at 240 to summer 242. Multiplier 244 multiplies 
the output 246 of C copy 146 and the error input at 44, 
and supplies the output resultant product at 248 to sum- 
mer 242. The summer 242 sums the inputs 248 and 240, 
and supplies the output resultant sum as weight update 
signal 78 to B filter 22. The input to A copy 214 and to 
B copy 216 is provided by the output 212 of stopband 
filter 21 0 receiving the noted random input signal at 206 
from random signal source 208. 

The error signals at error inputs 232 and 224 oppo- 
sitely drive the model. The error signal at error input 232 
of the direct transfer function filter A drives the correction 
signal at 46 towards a value matching the system input 
signal 6. The error signal at error input 224 of filter A 
drives the correction signal at 46 away from the noted 
matching value by driving the correction signal towards 
zero. As noted above, this is accomplished by using a 
copy 21 4 of the A filter and supplying the output of such 
copy as an error input to the adaptive model such that 
in attempting to drive the error input to zero, the model 
must drive its output to zero. The signal at error input 
224 is provided only in response to a given condition of 
a given parameter, e.g. when the frequency is outside 
the stopband of filter 210. The relative amplitudes of the 
input signals at error inputs 232 and 224 are adjusted 
such that the signal at error input 224 dominates when 
both are present, or the degree of dominance is adjusted 
to in turn adjust the amount of constrainment of perform- 
ance of the model so that the correction signal at 46 is 
driven towards zero but never reaches zero, such that 
there is still some modification and/or cancellation of the 
system input signal, though to a reduced degree. When 
the frequency is in the stopband of filter 210, there is no 



output at 212, and hence no input to A copy 214 and 
hence the latter's output is undefined, whereby error in- 
put 232 from error signal 44 from error transducer 16 
dominates and hence drives correction signal 46 to a 

5 value which matches the system input signal 6 to pro- 
vide modification and/or control of the latter. The error 
inputs 248 and 240 to the recursive transfer function fil- 
ter B of model 40 function comparably to error inputs 
232 and 224, respectively. Model 40 has a first error in- 
fo put provided at 232 and 248 from error transducer 16 
driving the output of the model towards a value matching 
the system input signal 6. Model 40 has a second error 
input at 224 and 240 selectively driving the output of 
model 40 away from such matching value and instead 

15 driving the correction signal 46 towards zero. 

Fig. 1 4 is similar to Fig. 1 3 and uses like reference 
numerals where appropriate to facilitate understanding. 
Summer 260 sums the output 230 of C copy 144 and 
the output 21 2 of stopband filter 21 0 which supplies the 

20 input to A copy 214, and supplies the output resultant 
sum at 262 to multiplier 264. Summer 266 sums the out- 
put 268 of A copy 214 and the error input at 44, and 
supplies the output resultant sum at 270 to multiplier 
264. Multiplier 264 multiplies the inputs 262 and 270, 

25 and supplies the output resultant product as weight up- 
date signal 74 to A filter 1 2. Summer 272 sums the out- 
put 246 of C copy 146 and the output 212 of stopband 
filter 210 which supplies the input to B copy 216, and 
supplies the output resultant sum at 274 to multiplier 

30 276. Summer 278 sums the output 280 of B copy 216 
and the error input at 44, and supplies the output result- 
ant sum at 282 to multiplier 276. Multiplier 276 multiplies 
inputs 282 and 274, and supplies the output resultant 
product as weight update signal 78 to B filter 22. 

35 The error signals at 44 and 268 oppositely drive the 
model. The error signal at error input 44 of the direct 
transfer function filter A drives the correction signal at 
46 towards a value matching the system input signal 6. 
The error signal at error input 268 of filter A drives the 

40 correction signal at 46 away from the noted matching 
value by driving the correction signal towards zero. As 
noted above, this is accomplished by using a copy 214 
of the A filter and supplying the output of such copy as 
an error input to the adaptive model such that in attempt - 

45 jng to drive the error input to zero, the model must drive 
its output to zero. The signal at error input 268 is pro- 
vided only in response to a given condition of a given 
parameter, e.g., when the frequency is outside the stop- 
band of filter 210. The relative amplitudes of the input 

50 signals at error inputs 44 and 268 to summer 266 are 
adjusted such that the signal at error input 268 domi- 
nates when both are present, or the degree of domi- 
nance is adjusted to in turn adjust the amount of con- 
strainment of performance of the model so that the cor- 

55 rection signal at 46 is driven towards zero but never 
reaches zero, such that there is still some modification 
and/or cancellation of the system input signal, though to 
a reduced degree. When the frequency is in the stop- 
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band of filter 210, there is no output at 212, and hence 
no input to A copy 214 and hence the latter's output is 
undefined, whereby error input 44 to summer 266 from 
error transducer 16 dominates and hence the model 
drives correction signal 46 to a value which matches the 
system input signal 6 to provide modification and/or con- 
trol of the latter. The error inputs 44 and 280 at summer 
278 to the recursive transfer function filter B of model 
40 function comparably to error inputs 44 and 268 at 
summer 266, respectively. Model 40 in Fig. 14 has a first 
error input provided at 44 to summers 266 and 278 from 
error transducer 16 driving the output of the model to- 
wards a value match ing the system input signal 6. Model 
40 has a second error input at 268 and 280 selectively 
driving the output of model 40 away from such matching 
value and instead driving the correction signal 46 to- 
wards zero. 

Fig. 15 illustrates the present invention and uses 
like reference numerals from above and from Figs. 19 
and 20 of the incorporated '676 patent. Random noise 
signal 206 from random noise source 208 is supplied 
through a first inverse C model copy 302 whose output 
is supplied to the filter input 304 of M copy 200. Correc- 
tion signal 46 from the output of model 40 through sum- 
mer 152 is supplied through a second inverse C model 
copy 306 to the convergence rate gain control input 308 
of M copy 200. The spectral transfer function provided 
by filter 302 spectrally controls performance of model 
40 to maximize the correction signal at frequencies of 
interest. In preferred form, the spectral transfer function 
is the inverse C model copy having the characteristics 
shown in Fig. 8, though other frequency responsive 
spectral transfer functions may be used. The spectral 
leak signal provided through transfer functions 302 and 
200 degrades performance of model 40, as above de- 
scribed in conjunction with Fig. 1 2. Transfer function 302 
controls the leak signal according to frequency. It is pre- 
ferred that the leak signal also be controlled in response 
to the correction signal from the output of model 40. 
Transfer function 306 spectrally monitors the correction 
signal and filters same to provide a spectral control sig- 
nal controlling the gain through M copy 200. Transfer 
function 306 and the gain control provided at control in- 
put 308 are optional. The spectral leak control provided 
through transfer function 302 to filter input 304 may be 
used with or without transfer function 306 and control 
input 308. Filter input 304 of M copy 200 receives a 
spectrally processed input signal from spectral transfer 
function 302. Control input 308 of M copy 200 controls 
gain of M copy 200 and receives a spectral control signal 
from transfer function 306. The spectral control signal is 
responsive to correction signal 46. Spectral transfer 
functions 302 and 306 may be the same, as shown, or 
may be different. 

Fig. 16 uses like reference numerals from above 
and from Figs. 1 9 and 20 of the incorporated '676 patent. 
As noted in the incorporated '676 patent, model M at 40 
is preferably an adaptive recursive filter having a trans- 



fer function with both poles and zeros. Model M is pro- 
vided by an MR, infinite impulse response, filter, e.g. a 
recursive least mean square, RLMS, filter having a first 
algorithm filter provided by an FIR, finite impulse re- 

5 sponse, filter, e.g. a least mean square, LMS, filter A at 
12, and a second algorithm filter provided by an FIR fil- 
ter, e.g. an LMS filter, B at 22. Filter A provides a direct 
transfer function, and filter B provides a recursive trans- 
fer function. The transfer function from output transduc- 

10 er 14 to error transducer 16 is modeled by a filter, e.g. 
an LMS or RLMS filter, C at 142, as in the incorporated 
'676 patent. 

Auxiliary random signal source 140 introduces a 
random signal into the output of model 40 at summer 

15 152 and into the C model at 148. The auxiliary random 
signal from source 140 is random and uncorrelated with 
the system input signal 6, and is also uncorrelated with 
auxiliary random signal source 208, and in preferred 
form is provided by a Galois sequence, M.R. Schroeder, 

20 "Number Theory In Science And Communications", Ber- 
lin, Springer-Berlag, 1 984, pages 252-261 , though other 
random uncorrelated signal sources may be used. The 
Galois sequence is a pseudo random sequence that re- 
peats after 2 M -1 points, where M is the number of stages 

25 in a shift register. The Galois sequence is preferred be- 
cause it is easy to calculate and can easily have a period 
much longer than the response time of the system. The 
input 1 48 to C model 1 42 is multiplied with the error sig- 
nal from error transducer 1 6 at multiplier 68, and the re- 

30 sultant product provided as weight update signal 67. 
Model 142 models the transfer function from output 
transducer 14to error transducer 16, including the trans- 
fer function of each. Alternatively, the transfer function 
from output transducer 14 to error transducer 16 may 

35 be modeled without a random signal source, as in U.S. 
Patent 4,987,598, incorporated herein by reference. 
Auxiliary source 1 40 introduces an auxiliary random sig- 
nal such that error transducer 1 6 also senses the auxil- 
iary signal from the auxiliary source. The auxiliary signal 

40 may be introduced into the recursive loop of the A and 
B filters as in Fig. 19 of the incorporated '676 patent at 
summer 1 52, or alternatively the auxiliary signal may be 
introduced into the model after the recursive loop, i.e. 
introducing the auxiliary signal only to line 46, and not 

45 to line 47. As in the incorporated '676 patent, copies of 
model 142 are provided at 144 and 146 to compensate 
the noted transfer function. 

The outputs of filters A and B are summed at sum- 
mer 48, whose output is summed at summer 152 with 

50 the output of random signal source 140 to provide an 
output resultant sum which provides the model output 
at 46 supplying the noted correction signal to output 
transducer 14. The output of model 142 is summed at 
summer 64 with the output of error transducer 16, and 

55 the resultant sum supplied as the error input to model 
142 and as an error input to model 40. Alternatively, the 
output of error transducer 16 may be supplied directly 
to an error input of model 40 without being supplied 
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through summer 64. 

M copy 200, Fig. 1 5, is provided by a copy of A filter 
12 at A copy 214, Fig. 16, and a copy of B filter 22 at B 
copy 216. Multiplier 218 multiplies the output 220 of A 
copy 21 4 and the input 222 of A copy 214, and supplies 
the output resultant product at 224 to summer 226. Mul- 
tiplier 228 multiplies the output 230 of C copy 144 and 
the error input at 44, and supplies the output resultant 
product at 232 to summer 226. Summer 226 sums the 
inputs 232 and 224, and supplies the output resultant 
sum as weight update signal 74 to A filter 12. Multiplier 
234 multiplies the output 236 of B copy 21 6 and the input 
238 of B copy 216, and supplies the output resultant 
product at 240 to summer 242. Multiplier 244 multiplies 
the output 246 of C copy 146 and the error input at 44, 
and supplies the output resultant product at 248 to sum- 
mer 242. The summer 242 sums the inputs 248 and 240, 
and supplies the output resultant sum as weight update 
signal 78 to B filter 22. The input to A copy 214 and to 
B copy 216 is provided by the output 304 of filter 302 
receiving the noted random input signal at 206 from ran- 
dom signal source 208. 

The error signals at error inputs 232 and 224 oppo- 
sitely drive the model. The error signal at error input 232 
of the direct transfer function filter A drives the correction 
signal at 46 towards a value matching the system input 
signal 6. The error signal at error input 224 of filter A 
drives the correction signal at 46 away from the noted 
matching value by driving the correction signal towards 
zero. As noted above, this is accomplished by using a 
copy 21 4 of the A filter and supplying the output of such 
copy as an error input to the adaptive model such that 
in attempting to drive the error input to zero, the model 
must drive its output to zero. The signal at error input 
224 is provided only in response to a given condition of 
a given parameter, e.g. when the frequency is in the 
ranges passed by filter 302. The relative amplitudes of 
the input signals at error inputs 232 and 224 are adjust- 
ed such that the signal at error input 224 dominates 
when both are present, or the degree of dominance is 
adjusted to in turn adjust the amount of constrainment 
of performance of the model so that the correction signal 
at 46 is driven towards zero but never reaches zero, 
such that there is still some modification and/or cancel- 
lation of the system input signal, though to a reduced 
degree. When the frequency is outside of the ranges 
passed by filter 302, there is no output at 21 2, and hence 
no input to A copy 214 and hence the latter's output is 
undefined, whereby error input 232 from error signal 44 
from error transducer 16 dominates and hence drives 
correction signal 46 to a value which matches the sys- 
tem input signal 6 to provide modification and/or control 
of the latter. The gain is controlled by the control signal 
at control input 308 which is the spectrally processed 
correction signal supplied through the transfer function 
filter provided by inverse C model copy 306. The error 
inputs 248 and 240 to the recursive transfer function fil- 
ter B of model 40 function comparably to error inputs 



232 and 224, respectively. Model 40 has a first error in- 
put provided at 232 and 248 from error transducer 16 
driving the output of the model towards a value matching 
the system input signal 6. Model 40 has a second error 
5 input at 224 and 240 selectively driving the output of 
model 40 away from such matching value and instead 
driving the correction signal 46 towards zero. 

Fig. 17 uses like reference numerals from above 
and from Figs. 1 9 and 20 of the incorporated '676 patent. 
Summer 260 sums the output 230 of C copy 144 and 
the output 21 2 of stopband filter 210 which supplies the 
input to A copy 214, and supplies the output resultant 
sum at 262 to multiplier 264. Summer 266 sums the out- 
put 268 of A copy 214 and the error input at 44, and 
supplies the output resultant sum at 270 to multiplier 
264. Multiplier 264 multiplies the inputs 262 and 270, 
and supplies the output resultant product as weight up- 
date signal 74 to A filter 1 2. Summer 272 sums the out- 
put 246 of C copy 146 and the output 304 of filter 302 
which supplies the input to B copy 21 6, and supplies the 
output resultant sum at 274 to multiplier 276. Summer 
278 sums the output 280 of B copy 216 and the error 
input at 44, and supplies the output resultant sum at 282 
to multiplier 276. Multiplier 276 multiplies inputs 282 and 
274, and supplies the output resultant product as weight 
update signal 78 to B filter 22. 

The error signals at 44 and 268 oppositely drive the 
model. The error signal at error input 44 of the direct 
transfer function filter A drives the correction signal at 
46 towards a value matching the system input signal 6. 
The error signal at error input 268 of filter A drives the 
correction signal at 46 away from the noted matching 
value by driving the correction signal towards zero. As 
noted above, this is accomplished by using a copy 214 
of the A filter and supplying the output of such copy as 
an error input to the adaptive model such that in attempt- 
ing to drive the error input to zero, the model must drive 
its output to zero. The signal at error input 268 is pro- 
vided only in response to a given condition of a given 
parameter, e.g., when the frequency is in the ranges 
passed by transfer function filter 302. The relative am- 
plitudes of the input signals at error inputs 44 and 268 
to summer 266 are adjusted such that the signal at error 
input 268 dominates when both are present, or the de- 
gree of dominance is adjusted to in turn adjust the 
amount of constrainment of performance of the model 
so that the correction signal at 46 is driven towards zero 
but never reaches zero, such that there is still some 
modification and/or cancellation of the system input sig- 
nal, though to a reduced degree. When the frequency 
is outside of the ranges passed by filter 302, there is no 
output at 212, there is no output at 212, and hence no 
input to A copy 214 and hence the latter's output is un- 
defined, whereby error input 44 to summer 266 from er- 
ror transducer 16 dominates and hence the model 
drives correction signal 46 to a value which matches the 
system input signal 6 to provide modification and/or con- 
trol of the latter. The gain is controlled by the control sig- 
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nal at control input 308 provided by the spectrally proc- 
essed correction signal supplied through the transfer 
function filter provided by inverse C model copy 306. 
The error inputs 44 and 280 at summer 278 to the re- 
cursive transfer function filter B of model 40 function 
comparably to error inputs 44 and 268 at summer 266, 
respectively. Model 40 in Fig. 6 has a first error input 
provided at 44 to summers 266 and 278 from error trans- 
ducer 16 driving the output of the model towards a value 
matching the system input signal 6. Model 40 has a sec- 
ond error input at 268 and 280 selectively driving the 
output of model 40 away from such matching value and 
instead driving the correction signal 46 towards zero. 

Fig. 18 illustrates the present invention and uses 
like reference numerals from above and from Figs. 19 
and 20 of the incorporated '676 patent. In Fig. 1 8, spec- 
tral transfer function filter 302 has its input 31 0 supplied 
from the correction signal from the output of model 40, 
rather than from random noise signal source 208, Fig. 
15. The spectral transfer function provided by filter 306 
also includes a peak detector 31 2 similar to peak detec- 
tor 97, Fig. 7. The output of inverse C model 306 is sup- 
plied to peak detector 312 comparing the output of in- 
verse C model 306 with a desired peak value 314. When 
the output of inverse C model copy 306 rises above the 
level of peak value 314, the positive output of summer 
312 at the control input 308 controls the gain of M copy 
200 to in turn control the gain of the leak signal at filter 
input 304 supplied through the spectral transfer function 
302. 

Fig. 19 uses like reference numerals from above 
and from Figs. 1 9 and 20 of the incorporated '676 patent. 
As noted in the incorporated '676 patent, model M at 40 
is preferably an adaptive recursive filter having a trans- 
fer function with both poles and zeros. Model M is pro- 
vided by an MR, infinite impulse response, filter, e.g. a 
recursive least mean square, RLMS, filter having a first 
algorithm filter provided by an FIR, finite impulse re- 
sponse, filter, e.g. a least mean square, LMS, filter A at 
12, and a second algorithm filter provided by an FIR fil- 
ter, e.g. an LMS filter, B at 22. Filter A provides a direct 
transfer function, and filter B provides a recursive trans- 
fer function. The transfer function from output transduc- 
er 14 to error transducer 16 is modeled by a filter, e.g. 
an LMS or RLMS filter, C at 142, as in the incorporated 
'676 patent. 

Auxiliary random signal source 140 introduces a 
random signal into the output of model 40 at summer 
152 and into the C model at 148. The auxiliary random 
signal from source 140 is random and uncorrelated with 
the system input signal 6, and is also uncorrelated with 
auxiliary random signal source 208, and in preferred 
form is provided by a Galois sequence, M.R. Schroeder, 
"Number Theory In Science And Communications"; Ber- 
lin, Springer-Berlag, 1 984, pages 252-261 , though other 
random uncorrelated signal sources may be used. The 
Galois sequence is a pseudo random sequence that re- 
peats after 2 M -1 points, where M is the number of stages 



in a shift register. The Galois sequence is preferred be- 
cause it is easy to calculate and can easily have a period 
much longer than the response time of the system. The 
input 1 48 to C model 1 42 is multiplied with the error sig- 
5 nal from error transducer 1 6 at multiplier 68, and the re- 
sultant product provided as weight update signal 67. 
Model 142 models the transfer function from output 
transducer 14to error transducer 16, including the trans- 
fer function of each. Alternatively, the transfer function 
10 from output transducer 14 to error transducer 16 may 
be modeled without a random signal source, as in U.S. 
Patent 4,987,598, incorporated herein by reference. 
Auxiliary source 1 40 introduces an auxiliary random sig- 
nal such that error transducer 1 6 also senses the auxil- 
15 jary signal from the auxiliary source. The auxiliary signal 
may be introduced into the recursive loop of the A and 
B filters as in Fig. 19 of the incorporated '676 patent at 
summer 1 52, or alternatively the auxiliary signal may be 
introduced into the model after the recursive loop, i.e. 
20 introducing the auxiliary signal only to line 46, and not 
to line 47. As in the incorporated '676 patent, copies of 
model 142 are provided at 144 and 146 to compensate 
the noted transfer function. 

The outputs of filters A and B are summed at sum- 
25 mer 48, whose output is summed at summer 152 with 
the output of random signal source 140 to provide an 
output resultant sum which provides the model output 
at 46 supplying the noted correction signal to output 
transducer 14. The output of model 142 is summed at 
30 summer 64 with the output of error transducer 16, and 
the resultant sum supplied as the error input to model 
142 and as an error input to model 40. Alternatively, the 
output of error transducer 16 may be supplied directly 
to an error input of model 40 without being supplied 
35 through summer 64. 

M copy 200, Fig. 1 8, is provided by a copy of A filter 
1 2 at A copy 21 4, Fig. 1 9, and a copy of B filter 22 at B 
copy 216. Multiplier 218 multiplies the output 220 of A 
copy 214 and the input 222 of A copy 214, and supplies 
40 the output resultant product at 224 to summer 226. Mul- 
tiplier 228 multiplies the output 230 of C copy 144 and 
the error input at 44, and supplies the output resultant 
product at 232 to summer 226. Summer 226 sums the 
inputs 232 and 224, and supplies the output resultant 
45 sum as weight update signal 74 to A filter 1 2. Multiplier 
234 multiplies the output 236 of B copy 21 6 and the input 
238 of B copy 216, and supplies the output resultant 
product at 240 to summer 242. Multiplier 244 multiplies 
the output 246 of C copy 1 46 and the error input at 44, 
50 and supplies the output resultant product at 248 to sum- 
mer 242. The summer 242 sums the inputs 248 and 240, 
and supplies the output resultant sum as weight update 
signal 78 to B filter 22. The input to A copy 214 and to 
B copy 216 is provided by the output 304 of filter 302 
55 receiving the model output at 47. 

The error signals at error inputs 232 and 224 oppo- 
sitely drive the model. The error signal at error input 232 
of the direct transfer function filter A drives the correction 
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signal at 46 towards a value matching the system input 
signal 6. The error signal at error input 224 of filter A 
drives the correction signal at 46 away from the noted 
matching value by driving the correction signal towards 
zero. As noted above, this is accomplished by using a 
copy 21 4 of the A filter and supplying the output of such 
copy as an error input to the adaptive model such that 
in attempting to drive the error input to zero, the model 
must drive its output to zero. The signal at error input 
224 is provided only in response to a given condition of 
a given parameter, e.g. when the frequency is in the 
ranges passed by filter 302. The relative amplitudes of 
the input signals at error inputs 232 and 224 are adjust- 
ed such that the signal at error input 224 dominates 
when both are present, or the degree of dominance is 
adjusted to in turn adjust the amount of constrainment 
of performance of the model so that the correction signal 
at 46 is driven towards zero but never reaches zero, 
such that there is still some modification and/or cancel- 
lation of the system input signal, though to a reduced 
degree. When the frequency is outside of the ranges 
passed by filter 302, there is no output at 21 2, and hence 
no input to A copy 214 and hence the latter's output is 
undefined, whereby error input 232 from error signal 44 
from error transducer 16 dominates and hence drives 
correction signal 46 to a value which matches the sys- 
tem input signal 6 to provide modification and/or control 
of the latter. The gain is controlled at control input 308 
by the control signal supplied from peak detector 312 
and filter 306 which supplies the spectrally processed 
correction signal based on the model output at 47. The 
error inputs 248 and 240 to the recursive transfer func- 
tion filter B of model 40 function comparably to error in- 
puts 232 and 224, respectively. Model 40 has a first error 
input provided at 232 and 248 from error transducer 16 
driving the output of the model towards a value matching 
the system input signal 6. Model 40 has a second error 
input at 224 and 240 selectively driving the output of 
model 40 away from such matching value and instead 
driving the correction signal 46 towards zero. 

Fig. 20 uses like reference numerals from above 
and from Figs. 1 9 and 20 of the incorporated '676 patent. 
Summer 260 sums the output 230 of C copy 144 and 
the output of filter 302 which supplies the input to Acopy 
21 4, and supplies the output resultant sum at 262 to mul- 
tiplier 264. Summer 266 sums the output 268 of A copy 
214 and the error input at 44, and supplies the output 
resultant sum at 270 to multiplier 264. Multiplier 264 mul- 
tiplies the inputs 262 and 270, and supplies the output 
resultant product as weight update signal 74 to A filter 
12. Summer 272 sums the output 246 of C copy 1 46 and 
the output 304 of filter 302 which supplies the input to B 
copy 216, and supplies the output resultant sum at 274 
to multiplier 276. Summer 278 sums the output 280 of 
B copy 216 and the error input at 44, and supplies the 
output resultant sum at 282 to multiplier 276. Multiplier 
276 multiplies inputs 282 and 274, and supplies the out- 
put resultant product as weight update signal 78 to B 



filter 22. 

The error signals at 44 and 268 oppositely drive the 
model. The error signal at error input 44 of the direct 
transfer function filter A drives the correction signal at 

s 46 towards a value matching the system input signal 6. 
The error signal at error input 268 of filter A drives the 
correction signal at 46 away from the noted matching 
value by driving the correction signal towards zero. As 
noted above ; this is accomplished by using a copy 214 

10 of the A filter and supplying the output of such copy as 
an error input to the adaptive model such that in attempt- 
ing to drive the error input to zero, the model must drive 
its output to zero. The signal at error input 268 is pro- 
vided only in response to a given condition of a given 

15 parameter, e.g., when the frequency is in the ranges 
passed by filter 302. The relative amplitudes of the input 
signals at error inputs 44 and 268 to summer 266 are 
adjusted such that the signal at error input 268 domi- 
nates when both are present, or the degree of domi- 
no nance is adjusted to in turn adjust the amount of con- 
strainment of performance of the model so that the cor- 
rection signal at 46 is driven towards zero but never 
reaches zero, such that there is still some modification 
and/or cancellation of the system input signal, though to 

25 a reduced degree. When the frequency is outside of the 
ranges passed by filter 302, there is no output at 212, 
and hence no input to A copy 21 4 and hence the latter's 
output is undefined, whereby error input 44 to summer 
266 from error transducer 1 6 dominates and hence the 

30 model drives correction signal 46 to a value which 
matches the system input signal 6 to provide modifica- 
tion and/or control of the latter. The gain is controlled at 
control input 308 by the control signal provided from 
peak detector 31 2 and filter 306. The error inputs 44 and 

35 280 at summer 278 to the recursive transfer function fil- 
ter B of model 40 function comparably to error inputs 44 
and 268 at summer 266, respectively. Model 40 in Fig. 
6 has a first error input provided at 44 to summers 266 
and 278 from error transducer 16 driving the output of 

40 the model towards a value matching the system input 
signal 6. Model 40 has a second error input at 268 and 
280 selectively driving the output of model 40 away from 
such matching value and instead driving the correction 
signal 46 towards zero. 

45 Fig. 21 uses like reference numerals from Fig. 18 
and illustrates a further embodiment. In Fig. 21 , spectral 
transfer function filter 302 has its input 31 6 supplied from 
the reference signal from input transducer 10, rather 
than from the correction signal from the output of model 

50 40, Fig. 18, or the random noise signal source 208, Fig. 
15. The output of spectral transfer function filter 302 is 
supplied to filter input 304 of M copy 200, as above. The 
spectral transfer function provided by filter 306 includes 
a peak detector 31 2 similar to peak detector 97, Fig. 7. 

55 The output of inverse C model 306 is supplied to peak 
detector 312 comparing the output of inverse C model 
306 with a desired peak value 314. When the output of 
inverse C model copy 306 rises above the level of peak 
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value 31 4, the positive output of summer 31 2 at the con- 
trol input 308 controls the gain of M copy 200 to in turn 
control the gain of the leak signal at filter input 304 sup- 
plied through the spectral transfer function 302 from the 
reference signal from input transducer 10. 

Filters 95, 302 ; 306 are each preferably an inverse 
C model copy provided by inverse S and/or inverse E, 
though other transfer functions may be used for any or 
all of such filters. 

In further embodiments, performance of model 40 
is controlled according to fuzzy logic to control the signal 
sent to output transducer 14. Fuzzy logic control is 
known in the prior art, for example "Adaptive Fuzzy Sys- 
tems", E. Cox, IEEE Spectrum February 1993, pages 
27-31 , and the references noted therein at page 31 , low- 
er half of right column. In the present invention, fuzzy 
logic is used to provide self-designing control architec- 
ture using fuzzy rules and/or to control the filter weights 
which update the model and/or to control the leak signal 
which degrades performance of the model. The fuzzy 
logic controller, having a given set of rules, is used for 
computing or setting the control architecture and/or the 
filtertransfer function and/or the filter weights and/or the 
leak signal. 

Fuzzy logic enables control of model performance 
in a practical way for a multiplicity of factors which inter- 
act in a complex way. The fuzzy logic control is based 
upon relative values and qualitative trends, without re- 
quiring exact equation relationships for multiple varia- 
bles. For example, referring to Figs. 22-24, which are 
similar to Fig. 4, page 29 of the above noted Cox article, 
if the magnitude of the weights is small, Fig. 22, and the 
magnitude of the correction signal output of model 40 is 
smalL Fig. 23, then only a small leakage value is intro- 
duced. Fig. 24, to maximize the adaptive process. If the 
magnitude of the weights is large and the magnitude of 
the model output is large, then a large leak signal is in- 
troduced, to constrain model performance and minimize 
the adaptive process. In this example, the filter weights 
are monitored to provide a first input parameter, and the 
output of the model is monitored to provide a second 
input parameter. Fuzzy logic is applied to such first and 
second input parameters to provide an output parame- 
ter, Fig. 24, controlling the weights and/or leak and/or 
otherwise controlling performance of the model. The in- 
put parameters are fuzzified according to a fuzzy rule 
set, Figs. 22 and 23, a fuzzy leak is computed and then 
def uzzified using Fig. 24, to provide an output parameter 
to control performance of the model. 

In another example, if the correction signal magni- 
tude is nearing the capacity of output transducer 1 4 such 
as a loudspeaker, then leakage is increased. If the mag- 
nitude of the correction signal is not nearing the capacity 
of the loudspeaker, then leakage can be decreased. 

In another example, rate of change of the weights 
is monitored, and if the weights start increasing at a fast- 
er rate or at a rate above a given rate, then leakage is 
increased. If the weights are increasing rapidly and the 



magnitude of the correction signal is approaching the 
capacity of the loudspeaker, then leakage is increased 
very rapidly. If the weights start decreasing, then leak- 
age is decreased. 
$ Inputs to the fuzzy logic include magnitude of the 

filter update weight, magnitude of the model output cor- 
rection signal, rate of change of the filter weights, rate 
of change of the correction signal output by the model, 
magnitude of the reference signal 42 input to the model, 
10 rate of change of input signal 42, magnitude of the error 
signal from error transducer 16, rate of change of the 
error signal, spectral characteristics of the reference sig- 
nal and/or correction signal and/or error signal, capabil- 
ity of output transducer 14, temperature, flow rates, en- 
's vironmental variables, fan speed in a duct application, 
whether or not such fan is running, thermostat settings, 
desired speed of adaptation, desired algorithm stability, 
system plant information, source information, etc. Fuzzy 
logic is a means of introducing human-like intuition into 
20 the controller adjustment process. Some factors may 
override others. For example, if the rate of increase of 
the weights is too large, then the algorithm may become 
unstable, and hence it is desirable to increase the leak 
or decrease the model update gain. If the capacity of a 
25 loudspeaker 14 is rapidly being reached, then leakage 
should be increased. 

In further embodiments, the invention is applicable 
to multi-channel active acoustic attenuation systems, for 
example as shown in U.S. Patents 5,216,721 and 
30 5,216,722, incorporated herein by reference. 

It is recognized that various equivalents, alterna- 
tives and modifications are possible within the scope of 
the appended claims. 
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Claims 



1 . An active adaptive control method comprising intro- 
ducing a control signal from an output transducer to 

40 combine with a system input signal and yield a sys- 
tem output signal, sensing said system output sig- 
nal with an error transducer providing an error sig- 
nal, providing an adaptive filter model having a 
model input from a reference signal correlated to 

45 said system input signal, and an output outputting 
a correction signal to said output transducer to in- 
troduce said control signal, spectrally controlling 
performance of said model to maximize the signal 
sent to said output transducer at frequencies of in- 

50 terest. 

2. The method according to claim 1 comprising pro- 
viding a spectral leak signal degrading performance 
of said model, and controlling said leak signal ac- 

55 cording to frequency. 

3. The method according to claim 2 comprising mon- 
itoring said correction signal and controlling said 
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leak signal in response thereto. 

4. The method according to claim 3 comprising filter- 
ing said correction signal to provide said leak signal. 

5. The method according to claim 4 comprising spec- 
trally processing said correction signal by supplying 
said correction signal through a frequency respon- 
sive spectral transfer function. 

6. The method according to claim 5 comprising mod- 
eling the transfer function between said output 
transducer and said error transducer with a second 
model, and wherein said spectral transfer function 
is a function of said second model. 

7. The method according to claim 5 comprising sup- 
plying said correction signal from said model output 
through said frequency responsive spectral transfer 
function to an error input of said model. 

8. The method according to claim 3 wherein said mod- 
el output outputs said correction signal to said out- 
put transducer to introduce said control signal ac- 
cording to a weight update signal, and comprising 
combining said reference signal and said error sig- 
nal to provide said weight update signal, and con- 
trolling said leak signal by controlling said weight 
update signal in response to said correction signal. 

9. The method according to claim 2 comprising pro- 
viding a copy of said model, said copy having an 
output supplied to said error input of said model, 
said copy having a filter input receiving a spectrally 
processed input signal from a spectral transfer func- 
tion, said copy having a control input controlling 
convergence gain of said copy and receiving a 
spectral control signal. 

10. The method according toclaim 9 wherein said spec- 
tral control signal is responsive to said correction 
signal. 

11. The method according to claim 10 comprising pro- 
viding a second spectral transfer function, and sup- 
plying said correction signal through said second 
spectral transfer function to said control input of said 
copy to provide said spectral control signal. 

12. The method according to claim 11 wherein said first 
and second spectral transfer functions are the 
same. 

13. The method according to claim 11 comprising mod- 
eling the transfer function between said output 
transducer and said error transducer with a second 
model, and wherein each of said first and second 
spectral transfer functions is a function of said sec- 



ond transfer function. 

14. The method according to claim 1 1 wherein said first 
and second spectral transfer functions are different. 

5 

15. The method according to claim 2 comprising spec- 
trally processing said leak signal such that 

at frequencies where maximum power from 
10 said output transducer reaches said error trans- 

ducer, the correction signal supplied to said 
output transducer is maximized, and 
at frequencies where minimum power from said 
output transducer reaches said error transduc- 
es er, the correction signal supplied to said output 
transducer is minimized. 

16. The method according to claim 1 5 comprising mod- 
eling the transfer function between said output 

20 transducer and said error transducer with a second 
model, and spectrally processing said leak signal 
as a function of said second model. 

17. The method according to claim 16 wherein said 
25 function of said second model is the inverse of said 

second model. 

18. The method according to claim 17 wherein said first 
model is controlled by an update signal from said 

30 inverse of said second model such that 

at frequencies where maximum power from 
said output transducer reaches said error trans- 
ducer, said second model has a maximum 
35 transfer characteristic and said inverse of said 

second model has a minimum transfer charac- 
teristic minimizing leakage of said update sig- 
nal, to enable maximum output of said first 
model, and 

40 at frequencies where minimum power from said 

output transducer reaches said error transduc- 
er, said second model has a minimum transfer 
characteristic and said inverse of said second 
model has a maximum transfer characteristic 
45 maximizing leakage of said update signal, to 

minimize the output of said first model. 

19. The method according to claim 2 comprising filter- 
ing said reference signal to provide said leak signal. 

50 

20. The method according to claim 1 9 comprising spec- 
trally processing said reference signal by supplying 
said reference signal through a frequency respon- 
sive spectral transfer function. 

55 

21 . The method according to claim 20 comprising mod- 
eling the transfer function between said output 
transducer and said error transducer with a second 
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model, and wherein said spectral transfer function 
is a function of said second model. 

22. The method according to claim 1 9 comprising mon- 
itoring said correction signal and controlling said 
leak signal in response thereto. 

23. The method according to claim 1 comprising spec- 
trally processing said leak signal by supplying said 
leak signal through a frequency responsive spectral 
transfer function. 

24. The method according to claim 23 wherein said leak 
signal is provided by said correction signal, and 
comprising supplying said correction signal through 
said frequency responsive spectral transfer func- 
tion. 

25. The method according to claim 24 comprising sup- 
plying said correction signal from said model output 
through said frequency responsive spectral transfer 
function to said error input of said model. 

26. The method according to claim 25 comprising pro- 
viding a copy of said model, and supplying said cor- 
rection signal from said model output through said 
frequency responsive spectral transfer function and 
through said copy to said error input of said model. 

27. The method according to claim 23 comprising pro- 
viding a copy of said model, and supplying said leak 
signal through said frequency responsive spectral 
transfer function and through said copy to said error 
input of said model. 

28. The method according to claim 27 wherein said 
copy has a filter input receiving the output of said 
frequency responsive spectral transfer function, 
said copy has a control input controlling conver- 
gence gain of said copy and comprising providing 
a second frequency responsive spectral transfer 
function, and supplying said correction signal from 
said output of said model through said second fre- 
quency responsive spectral transfer function to said 
control input of said copy. 

29. The method according to claim 27 comprising sup- 
plying a random noise signal to the input of said 
spectral transfer function. 

30. The method according to claim 27 comprising sup- 
plying said correction signal from said output of said 
model to the input of said spectral transfer function. 

31. The method according to claim 1 comprising spec- 
trally controlling performance of said model in re- 
sponse to said correction signal. 



32. The method according to claim 31 wherein said 
model outputs said correction signal to said output 
transducer to introduce said control signal accord- 
ing to a weight update signal, and comprising com- 

s bining said reference signal and said error signal to 

provide said weight update signal, monitoring and 
spectrally sensing said correction signal and provid- 
ing selective leakage of said weight update signal 
in response thereto to control performance of said 

10 model according to frequency, to optimize perform- 
ance of said model in frequency ranges of interest. 

33. The method according to claim 31 comprising pro- 
viding a copy of said model and supplying the output 

15 of said copy to said error input of said model and 
monitoring and spectrally sensing said correction 
signal and providing an input to said copy in re- 
sponse thereto to control performance of said mod- 
el according to frequency to optimize performance 

20 of said model in frequency ranges of interest. 

34. The method according to claim 2 comprising fre- 
quency weighting said leak signal to optimize per- 
formance of said model in frequency ranges of in- 

25 terest. 

35. The method according to claim 34 comprising fre- 
quency weighting said leak signal by spectrally 
processing said correction signal through a fre- 

30 quency responsive spectral transfer function. 

36. The method according to claim 34 comprising fre- 
quency weighting said leak signal by spectrally 
processing said correction signal through a fre- 

35 quency responsive spectral transfer function and 
through a copy of said model. 

37. The method according to claim 34 comprising fre- 
quency weighting said leak signal by spectrally 

40 processing a random noise signal through a fre- 
quency responsive spectral transfer function and 
through a copy of said model. 

38. The method according to claim 34 comprising fre- 
45 quency weighting said leak signal by spectrally 

processing a random noise signal through a first fre- 
quency responsive spectral transfer function and 
through a copy of said model, and spectrally 
processing said correction signal through a second 
50 frequency responsive spectral transfer function and 
through said copy of said model. 

39. The method according to claim 34 comprising fre- 
quency weighting said leak signal by spectrally 

55 processing said correction signal through a first fre- 
quency responsive spectral transfer function and 
through a copy of said model, and spectrally 
processing said correction signal through a second 
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frequency responsive spectral transfer function and 
through said copy of said model. 

40. The method according to claim 1 comprising spec- 
trally controlling performance of said model in re- 
sponse to said reference signal. 

41. The method according to claim 40 wherein said 
model outputs said correction signal to said output 
transducer to introduce said control signal accord- 
ing to a weight update signal, and comprising com- 
bining said reference signal and said error signal to 
provide said weight update signal, monitoring and 
spectrally sensing said reference signal and provid- 
ing selective leakage of said weight update signal 
in response thereto to control performance of said 
model according to frequency, to optimize perform- 
ance of said model in frequency ranges of interest. 

42. An active adaptive control system comprising an 
output transducer introducing a control signal to 
combine with a system input signal and yield a sys- 
tem output signal, an error transducer sensing said 
system output signal and providing an error signal, 
an adaptive filter model having a model input from 
a reference signal correlated to said system input 
signal, and an output outputting a correction signal 
to said output transducer to introduce said control 
signal, a spectral controller spectrally controlling 
performance of said model to maximize the signal 
sent to said output transducer at frequencies of in- 
terest. 

43. An active adaptive control method comprising intro- 
ducing a control signal from an output transducer to 
combine with a system input signal and yield a sys- 
tem output signal, sensing said system output sig- 
nal with an error transducer providing an error sig- 
nal, providing an adaptive filter model having a 
model input from a reference signal correlated to 
said system input signal, and an output outputting 
a correction signal to said output transducer to in- 
troduce said control signal, controlling performance 
of said model according to fuzzy logic to control the 
signal sent to said output transducer. 

44. The method according to claim 43 comprising pro- 
viding a leak signal degrading performance of said 
model, and controlling said leak signal according to 
said fuzzy logic. 

45. The method according to claim 43 comprising up- 
dating said model by a set of filter weights, and con- 
trolling said weights according to said fuzzy logic. 

46. The method according to claim 43 comprising f uzz- 
ifying a plurality of input parameters according to a 
fuzzy rule set, computing a fuzzy leak, and def uzz- 



ifying said leak to provide an output parameter to 
control performance of said model. 

47. The method according to claim 43 comprising up- 
$ dating said model by a set of filter weights, monitor- 
ing said weights to provide a first input parameter, 
monitoring the output of said model to provide a 
second input parameter, applying fuzzy logic to said 
first and second input parameters to provide an out- 

10 put parameter controlling said weights. 

48. The method according to claim 47 comprising pro- 
viding a leak signal degrading performance of said 
model, and controlling said leak signal with said out- 

15 put parameter. 

49. An active adaptive control method comprising intro- 
ducing a control signal from an output transducer to 
combine with a system input signal and yield a sys- 

20 tern output signal, fuzzifying a plurality of input pa- 
rameters according to a fuzzy rule set and providing 
a control architecture therefrom outputting a correc- 
tion signal to said output transducer to introduce 
said control signal. 

25 

50. The method according to claim 49 comprising pro- 
viding said control architecture by setting a filter 
transfer function. 

30 51. The method according to claim 50 comprising set- 
ting said filter transfer function by setting filter 
weights. 



52. The method according to claim 51 comprising pro- 
35 viding a leak signal degrading performance of said 

filter transfer function, computing a fuzzy leak from 
said fuzzified input parameters and then defuzzify- 
ing said leak to provide an output parameter to con- 
trol performance of said filter transfer function. 

40 

53. An active adaptive control system comprising an 
output transducer introducing a control signal to 
combine with a system input signal and yield a sys- 
tem output signal, an error transducer sensing said 

45 system output signal and providing an error signal, 
an adaptive filter model having a model input from 
a reference signal correlated to said system input 
signal, and an output outputting a correction signal 
to said output transducer to introduce said control 

so signal, fuzzy logic controlling performance of said 
model to control the signal sent to said output trans- 
ducer. 

54. An active adaptive control method comprising intro- 
55 ducing a control signal from an output transducer to 

combine with a system input signal and yield a sys- 
tem output signal, sensing said system output sig- 
nal with an error transducer providing an error sig- 
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nal, providing an adaptive filter model having a 
model input from a reference signal correlated to 
said system input signal, and a model output out- 
putting a correction signal to said output transducer 
to introduce said control signal according to a 5 
weight update signal, adaptively varying said 
weight update signal. 

55. The method according to claim 54 comprising adap- 
tively varying said weight update signal as a func- 10 
tion of said correction signal and supplying the 
adaptively varied weight update signal to said adap- 
tive filter model to adapt said correction signal var- 
iably controlling leakage of said weight update sig- 
nal. 15 

56. The method according to claim 55 comprising var- 
ying leakage of said weight update signal by multi- 
plying a previous weight update value by a factor y 
and adding the result to the product of said refer- 20 
ence signal and said error signal, and varying y as 

a function of said correction signal. 

57. The method according to claim 54 comprising adap- 
tively varying said weight update signal as a func- 25 
tion of said correction signal relative to a given 
threshold to alter convergence of said correction 
signal when the latter exceeds said threshold. 
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